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MARK EMBEDDING AND DETECTION USING PROJECTIVE TRANSFORMS 

Field of the Invention 

4 

The present invention relates generally to the embedding and detection of 
5 watermarks in an image, and, in particular, to the embedding and detection of patterns having 
variation in one direction only. 

Background 

In the last decade numerous schemes have been proposed for the imperceptible 
embedding and detection of marks in signals. For a mark to be imperceptible in a signal, such 

10 as an image or even a video, the mark is typically embedded at a very low level when 
compared with the signal. For such a low level signal to' be detectable, a highly efficient 
detection process is needed. Such a process exists in the fonn of correlation or "matched 
filtering" in the case where the process is an optimal linear detection. 

Certain marks have near perfect correlation properties, allowing virtually all the 

15 energy of the mark to be concentrated into a single detection point. Unfortunately most 
correlation processes are effectively destroyed, or at least severely impeded, by distortions of 
the marked signal. In two dimensions, common distortions include rotation, scaling, 
anamorphic scaling and shear. 

Summary 

20 It is an object of the present invention to substantially overcome, or at least 

ameliorate, one or more disadvantages of existing arrangements. 

Disclosed are arrangements which seek to detect an embedded watermark, the 
watermark being formed from a one-dimensional function, even when the image including the 
watermark has been rotated, by applying a projective transform on the image, followed by 

25 correlation with the basis function. 
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According to an aspect of the present invention there is provided a method of 
detecting one or more patterns embedded in an image, each pattern having been fontied 
substantially from a one-dimensional basis function, said method comprising the steps of: 
(a) calculating a projective trans fomi of said image; 
5 (b) calculating a 1-D correlation betwreen the projective transform and said 

basis function for a selection of angles; and 

(c) finding one or more peaks of said correlation, wherein the position of each 
of said peaks provides spatial parameters of one of said one or more embedded patterns. 

Another aspect of the present invention provides a method of determining 
10 transformations applied to an image, said image having at least three patterns embedded 
therein, each pattern having been formed substantially from a one-dimensional basis function 
being substantially scale invariant as herein defined, said method comprising the steps of: 

(a) calculating a projective transform of said image; 

(b) calculating a 1-D correlation between the projective transform and said 
1 5 basis function for a selection of angles; 

(c) finding peaks of said correlation, wherein the positions of said peaks 
provide spatial parameters of said embedded patterns; and 

(d) determining fi-oni said spatial parameters said transformations. 

Another aspect of the present invention provides a method of detecting one or 
20 more patterns embedded in an image, each pattern having been formed substantially from a 
one-dimensional basis function, said method comprising the steps of: 

transforming said image to the frequency domain using a two-dimensional 

transform; 

resampling the transformed image onto a quasi-polar map; 
25 transforming said basis function to the frequency domain using a one-dimensional 

transform; 
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multiplying the resampled transformed image with the transform of said basis 
function along radial lines to create a first result; 

inverse transforming said first result lo the space domain using a one-dimensional 
transform to create a second result; and 
5 finding one or more peaks of said second result, wherein the position of each of 

said peaks provides spatial parameters of one of said one or more embedded patterns. 

Yet another aspect of the present invention provides a method of determining 
transformations applied to a first image, said first image being a transformed copy of a second 
image, said second image having at least three patterns embedded therein, each pattern having 
10 predefined parameters and having been fonned substantially from a one-dimensional basis 
function being substantially scale invariant as herein defined, said method comprising the 
steps of: 

calculating a projective transform of said first image; 

calculating a 1-D correlation between the projective transform and said basis 
15 function for a selection of angles; 

finding peaks of said correlation, each peak corresponding with one of said 
embedded patterns; and 

determining said transformations from the positions of said peaks and said 
parameters of said patterns embedded into said second image. 
20 Yet another aspect of the present invention provides a method of embedding a 

watermark into an image, said method comprising the steps of: 

maintaining a basis function, said basis function being a scale invariant one- 
dimensional function excluding the function /(r) = cos(aloglr| + c) wherein a and c are 

constants; 

25 forming one or more patterns from said basis function, each pattern having variation in 

one direction; and 
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adding said basis pattem(s) to said image. 

Oilier aspects of the invention are also disclosed. 

Brief Description of the Drawings 
One or more embodiments of the present invention will now be described wixh 
5 reference to the drawings, in which: 

Fig. lA shows a flow diagram of a method of detecting one or more patterns 
embedded in a two-dimensional function or signal; 

Fig, IB shows a flow diagram of a method, which is an implementation of the 
method shown in Fig. I A of detecting one or more patterns embedded in an image; 
10 Fig. IC shows a flow diagram of a method, which is another implementation of the 

method shown in Fig. 1 A; 

Fig. 2 shows a flow diagram of a method of embedding one or more patterns in the 

two-dimensional function; 

Fig. 3 shows a schematic block diagram of an apparatus for performing the 

1 5 operations of the above methods; 

Fig. 4A illustrates an example pattern; 

Fig. 4B illustrates the axis of symmetry of the pattern illustrated in Fig. 4A together 
with the orientation angle and perpendicular displacement from the centre of the pattern 

image; 

20 Fig. 4C illustrates a relationship between the vectors r and p; 

Fig. 5 illustrates the operation of the projective transform; 

Fig. 6 illustrates the projection pdr) of the example pattern illustrated in Fig. 4A 

plotted in the quasi polar space; 

Fig. 7 illustrates a resulting correlation of the 1-D correlation between the projection 
25 p6(r) illustrated in Fig. 6 and an associated basis function; 



583l88.doc 



-5- 



Fig. 8A illustrates 3 correlation peaks detected in an image having 3 embedded 
patterns with different parameters; 

Fig. 8B illustrates the correlation peaks detected from a distorted version of the same 

image as that in used for Fig. 8A; 
5 Fig. 9 shows a comparison between a fractional power-phase function ^ 

rescaled version ;i^(0.8r) thereof, as well as a phase mismatch between the functions; 

Fig. 10 shows a plot of how the magnitude of the cross-correlation of fractional 
power-phase ftinclions zi^) and ;i^(0.8r) as a function the fractional power and 

Fig. 1 1 shows a schematic block diagram of a system that may be used to implement 
10 the embedding and detection of watermarks in an image. 

Detailed Description including Best Mode 

Where reference is made in any one or more of the accompanying drawings to steps 
and/or features, which have the same reference numerals, those steps and/or features have for 
the purposes of this des.cription the same function(s) or operation(s), unless the contrary 

1 5 intention appears. 

Fig. 1 A shows a flow diagram of a method 200 of detecting one or more patterns 
embedded in a two-dimensional function or signal, such as an image. Before describing the 
method 200 of detecting the patterns, a method 300 of embedding such patterns in the two- 
dimensional function will be described with reference to Fig. 2. 

20 Some portions of the description which follows are explicitly or implicitly presented 

in terms of algorithms and symbolic representations of operations on data within a computer 
memory. These algorithmic descriptions and representations are the means used by those 
skilled in the data processing arts to most effectively convey the substance of their work to 
others skilled in the art. An algorithm is here, and generally, conceived to be a self-consistent 

25 sequence of steps leading to a desired result. The steps are those requiring physical 
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manipulations of physical quantities. Usually, though not necessarily, these quantities take 
the form of electrical or magnetic signals capable of being stored, transferred, combined, 
compared, and otherwise manipulated within registers and memories of a computer system, or 
similar electronic device. 

5 Fig. 3 shows a schematic block diagram of an apparatus in the form of a general- 

purpose computer system 100 for performing the operations of methods 200 and 300, wherein 
the steps of methods 200 and 300 may be implemented as software, such as one or more 
application program executing within the computer system 100. 

In particular, the steps of methods 200 and 300 are effected by instructions in the 

10 software that are carried out by the computer system 100. The software may be stored in a 
computer readable medium, including the storage devices described below, for example. The 
software is loaded into the computer system 100 from the computer readable medium, and 
then executed by the computer system 100. The use of the software in the computer system 
100 effects an advantageous apparatus for detecting one or more patterns embedded in a two- 

15 dimensional function, and an advantageous apparatus for embedding such patterns in the two- 
dimensional function. 

The computer system 100 is formed by a computer module 101, input devices 
including a keyboard 102, a mouse 103 and an imaging device 122, and output devices 
including a printer 1 15 and a display device 1 14. The imaging device 122 may be a scanner 

20 or digital camera used for obtaining a digital image. A Modulator-Demodulator (Modem) 
transceiver device 1 16 is used by the computer module 101 for communicating to and from a 
communications network 120, for example connectable via a telephone line 121 or other 
functional medium. The modem 1 1 6 can be used to obtain access to the Internet, and other 
network systems, such as a Local Area Network (LAN) or a Wide Area Network (WAN), and 

25 may be incorporated into the computer module 101 in some implementations. 
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The computer module 101 typically includes at least one processor unit 105, and a 
memory unit 106, for example fonned from semiconductor random access memory (RAM) 
and read only memory (ROM). The module 101 also includes a number of input/output (I/O) 
interfaces including a video interface 107 thai couples to the video display 114, an I/O 
5 interface 11 3 for the keyboard 102, mouse 103 and imaging device 122, and an interface 108 
for the modem 1 16 and printer 1 15. A storage device 109 is provided and typically includes a 
hard disk drive 110 and a floppy disk drive 111. A CD-ROM drive 1 1 2 is typically provided 
as a non-volatile source of data. The components 105 to 1 13 of the computer module 101, 
typically communicate via an interconnected bus 104 and in a manner which results in a 
10 conventional mode of operation of the computer system 100 known to those in the relevant 
art. 

Typically, the application program is resident on the hard disk drive 110 and read 
and controlled in its execution by the processor 105. Intermediate storage of the program and 
any data fetched from the network 120 may be accomplished using the semiconductor 

15 memory 106, possibly in concert with the hard disk drive 110. In some instances, the 
application program may be supplied to the user encoded on a CD-ROM or floppy disk and 
read via the corresponding drive 1 12 or 1 1 1, or alternatively may be read by the user from the 
network 120 via the modem device 1 16. Still further, the software can also be loaded into the 
computer system 100 from other computer readable media. The term "computer readable 

20 medium" as used herein refers to any storage or transmission medium that participates in 
providing instructions and/or data to the computer system 100 for execution and/or 
processing. Examples of storage media include floppy disks, magnetic tape, CD-ROM, a 
hard disk drive, a ROM or integrated circuit, a magneto-optical disk, or a computer readable 
card such as a PCMCIA card and the like, whether or not such devices are internal or external 

25 of the computer module 101. Examples of transmission media include radio or infra-red 
transmission channels as well as a network connection to another computer or networked 
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device, and the Internet or Intranets including e-mail transmissions and information recorded 
on Website? iand the hke. 

The methods 200 and 300 may alternatively be implemented in dedicated hardware 
such as one or more integrated circuits performing the functions or sub functions of methods 
5 200 and 300. Such dedicated hardware may include graphic processors, digital signal 
processors, or one or more microprocessors and associated memories. In one implementation, 
such integrated circuits are included in an imaging device, such as a camera or scanner. 

Referring again to Fig. 2, the method 300 of embedding patterns in the two- 
dimensional function starts in step 310 where a predetermined number of patterns are 
10 generated, each with predetermined parameters, and the separate patterns are summed 
together to form a pattern image. Such operations may be performed by the processor 105 
under control of one or more application programs. The pattern image may be stored after 
generation in the memory 106 or HDD 110. 

Each pattern has a variation in one direction only, that direction being at an 
15 orientation angle Pn with the horizontal direction. The pattern also has an axis of symmetry, 
which is perpendicular to the variation direction, and at a displacement v„ from the centre of 
the pattern image. The pattern may have even or odd symmetry around the axis of symmetry. 

The patterns are typically generated from a one-dimensional basis function applied at 
the direction of variation and repeated along the axis of symmetry. The basis function is 
20 preferably attenuated in areas where the basis function has a frequency above the Nyquist 
frequency of the pattern image. 

In a preferred implementation, the basis function is a complex homogeneous function 
of the form: 

g(r) = |r-v.„r" =|/--v„|"exp(mlog|/--vJ (1) 
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where r is a one-dimensional coordinate, which is not necessarily positive, Vq is an 
offset, while a and p are constants. Equation (1) may be considered as a function amplitude 
modulating a phase function, with the phase function having a logarithmic phase. When such 
a complex homogeneous function is scaled, say by a factor a, the scaling only introduces a 
5 complex constant factor as follows; 



The advantage of the complex homogeneous function is that the auto-correlation of 
the complex homogeneous function is directly proportional to the cross-correlation of the 
complex homogeneous function with a scaled version of the complex homogeneous function. 
10 This 'scale-invariant' properly allows a watemiark to be detected in an image even after a 
scale transformation by correlating the image with the basis function. 

In practice, a basis function useable as a scale invariant function only has to 
approximate this ideal property sufficiently so as to provide a dominant correlation peak. For 
example, a fractional power-phase function defined by: 



(2) 



!5 




(3) 



has a phase derivative (i.e. instantaneous frequency) that is similar to that of the ideal 



homogeneous function for small fractional powers in the rangc|y|<l, y^O. The phase 



derivative of the fractional power-phase function is: 



dr 



-or 



(4) 



20 



while the phase derivative of the ideal homogeneous function g(r) is: 



rfarg[g(r)] ^., 
dr 



(5) 
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Fig. 9 shows the real part of the fractional power-phase function ^(r) in the range 
(0.5. 1) as a continuous curve, and a version of the fractional power-phase function z{^') 
rescaled by a factor of 0,8 as a dotted curve. The fractional power-phase function xi'') has 
parameters p=0, a=606n, and j^O.l. The phase mismatch between the functions is also 
5 shown. The cross-correlation between the two functions illustrated will thus not be an ideal 
delta function, but rather some spread of energy around a point. 

Fig. 10 shows a plot of how the magnitude of the cross-correlation of fractional 
power-phase functions ^(r) and yi:(0.8r) varies as a function of the fractional power y. It can 
be observed from the plot that the correlation drops from about 70% at 7 = 0.1 to near 0% at 
10 y = Q.2. 

Generally, a measure of a function's scale invariance is equivalent to how much the 
phase profile of the function matches after scaling. 

The ideal cross-correlation or overlap integral for a nominally orthogonal 1-D 
function with inverse square root weighting, as is the case for the function in Equation (3) 
15 with /?=-0.5, is as follows: 



1 / / ^ f-^xp[/>(/0>/r 



(6) 



wherein i//{r) is a phase mismatch function. For two perfect phase matched 
functions the phase mismatch function v^(a-)=0, and the correlation is not attenuated. For 
small phase errors, that is the phase mismatch function \if/{xl < 1 , the cross-correlation energy 
20 is as follows: 



_L_'fi 



dx 



(7) 



Ignoring third and higher powers, and defining weighted phase statistics as: 
S83188.doc 



- II - 



ln(.v,/.v,)J A- 



(8) 



then Equation (7) may be simplified to: 



(9) 



Ignoring the last term, this may be expressed as: 



(10) 



wherein a is the weighted root mean square phase mismatch. It is noted that other 
weighting, that is other than l/Vr , may be applied to the phase statistics in Equation (8) if the 
need arises. 

Preferably, the phase profiles of the unsealed and the scaled basis function match 
10 (apart from a constant or linear tenn) within a RMS of 0.7 radian. In such a case, the cross- 
conrelation will not be attenuated by more than 30% when compared with the ideal. Table 1 
sets out the magnitude of the correlation for different values of the RMS phase mismatch a , 
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Table 1 



RMS phase 
misinatch cr 


Correlation 

|c„|a= Vl-O-' 


n 


1 .uu 


0 1 

\j, I 


n 00 
u, w 






U.J 




0.4 


0.92 


0.5 


0.87 


0.6 


0.80 


0.7 


0.71 


0.8 


0.60 


0.9 


0.44 



Other basis functions include: 

/, (/•) = [sgn(/-)]* exp(f (or/--" +£:)); (11) 

5 /2('-)=i;^,/oUv);aiid (12) 

«=/V, 

f, (r) = [sgn(r)]V'' expj^z ^r^ j (13) 

where e^p.Nx^ Niy ^. * and // are constants and the function is an arbitrary 
function. The constant fu is preferably in the range 0.5 to 2, but excluding 1. The constant 
k=Q or 1. In the case where ^=0, the function 7^(0 has even symmetry, whereas in the case 
10 where the function /^(r) has odd symmetry. 

Fig. 4A illustrates an example pattern generated using a real part of the basis function 
of Equation (1), with masking applied at the area adjoining the axis of symmetry to remove 
values with a frequency above the Nyquist frequency of the pattern image. In illustrating the 
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example pattern, pattern values has been mapped to values in the range of [0, 255], with a 
value of 0 being represented as the colour black and a value of 255 being represented as the 
colour while. Fig. 4B illustrates the axis of symmetry of the pattern illustrated in Fig. 4A 
together with the orientation angle (in and peipendicular displacement v,^ from the centre of 
5 the pattern image. 

Referring again to Fig. 2, in ord^r to embed patterns that are imperceptible to a 
human observer, the pattern image is retrieved from the memory 106 and is perceptually 
masked in step 320 by the processor 105 in order to greatly reduce the levels of the patterns 
corresponding to regions of the image having low intensity variation, and reduce by a lesser 

10 amount the levels of the patterns corresponding to regions of the image having high intensity 
variance. An example measure of intensity variation is the local gradient magnitude of the 
luminance in the image. Other measures include second partial derivatives of the luminance; 
local estimates of the "energy" or frequency content, local variance, and more sophisticated 
estimates of human visual system masking. 

15 The perceptually masked pattern image, which may be called a watermark, is added 

to the two-dimensional function, such as an image, in step 330. If the image is a colour 
image, then the watermark is preferably added to the luminance part of a colour image. This 
allows the watermark to survive when the watermarked image is converted from colour to a 
greyscale representation. Alternatively, the watermark may be added to one of the R, G, B, 

20 H, V, S, u, v etc channels of the colour image, or any combination thereof. Apart from simple 
algebraic a<idition, addition of the watermark to the image also includes dithering and half- 
toning. Watermarks may also be applied to sequences of images such as video streams by 
sequential addition. 

The image with the embedded watermark may be stored on storage device 109 (Fig. 
25 3) in a digital image format, such as tiff, gif, jpeg, mpeg etc. The watermarked image may 

583188.doc 



- 14- 



also be printed on printer 1 15 (Fig. 3) to create a hardcopy of the watermarked image, or may 
be communicated to another computer module 101 using the network 120. 

Referring again to Fig. lA where the flow diagram of the method 200 of detecting 
one or more patterns embedded in an image is shown, with the patterns embedded in the 
image using the method 300 (Fig. 2) described above. Method 200 starts in step 202 where a 
digital copy of the image is obtained. Typically the image may already be in a digital image 
format and stored on storage device 109 (Fig. 3). In such cases the image is converted to an 
array of pixel values. If the image is still in a hardcopy format, such as a photograph, then the 
hardcopy image is first converted into a digital format by scanning the hardcopy using the 
imaging device 122 (Fig. 3). 

Step 205 follows where the processor 105 undoes the perceptual masking by first 
forming a perceptual mask from the image, and then emphasising the image with the 
perceptual mask by dividing the values of the image by the corresponding values of the 
perceptual mask. 

A projective transform is then applied to the resulting image in step 210. The 
projective transform accumulates energy by summing values along straight lines in the image. 
This is similar to the way that the total effective path of X-rays accrues during transmission 
through a patient's body. Fig. 5 illustrates the operation of the projective transform in that a 
projection of the values of a function h(x^) is accumulated along straight lines with an angle 
0 with the vertical direction to form a projection p^r). The projection p^r) is found for a 
selection of projection angles ^in the range (-;r,;T]. 

The Radon (or equivalently Hough) transform is one such projective transform that 
may be used in step 210 and is defined as: 

^9 y)] = Ptf ('■) = J //»( V, y)S{r -xcosa-y sin 0)dxdy ( 1 4) 
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In order to derive a convenient implementation of the Radon transform for a discrete 
dataset, a useful correspondence between the projection p^r) of 2'D function h{x^f) and the 
slices of the function's Fourier transform is used, that correspondence being known as the 
"projection-slice theorem". 

The projection-slice theorem states that the one-dimensional Fourier Iransfomi of a 
projection of a two dimensional function is equal to a radial slice of the two-dimensional 
Fourier transform of that function. Hence: 

^(".v)= J \h{x,y)Q\p[-2m{nx + v\>y^dy (15) 

-09-00 

wherein H{u,v) is the 2-D Fourier transform of function h{x,y). In the quasi polar 
space the angles are in the range Tt/l^ir/l], while distance is in the range (-00,00). By 
defining quasi-polar coordinates {q,(p) in the Fourier domain, the coordinate transform is 

// = 7cos^, v = ^/sinj*, and one form of the projection-slice theorem is obtained for the 
Fourier polar angle corresponding to the Radon projection angle <^ = ^ as follows: 
+00 

\p0{r)e^^-2m{rq)\lr = H{q cos q s\ne)= Pg{q) (16) 

-GO 

Equation (16) is useful because it allows estimation of (the Fourier transform of) a 
Radon projection as a radial slice of the 2-D-FFT of a discrete image. This suggests that a 
discrete Radon transform may be evaluated by first performing a 2-D FFT followed by a 
Cartesian to polar remapping, using a suitable interpolation - such as bicubic, or chirp-z - to 
perform the resampling. L 

A two-dimensional function h{x^) formed from an one-dimensional basis function 
g(r) having an orientation angle fi„ and perpendicular displacement Vn may be written as: 

h{x,y) = g{xcosP^ +3^sin>3,, -v/J (17) 
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The two-dimensional function /»(a\)) formed from the complex homogeneous 
function of Equation (1) may be written in an even fonn and an odd form as follows: 

/w(^.3^) = |-vcos A, + J^sin A, -v/J^"" (18) 
/w(^'>^)=sgnCvcos>ff„ +.vsin/?„ -v/Ja-cos/?„ +.vsin^„ -^T"" (19) 

Linear combinations of the functions in Equations (18) and (19) may also be used. 
Note that a one-sided function is obtained from adding an odd and even function in equal 
proportion. 

When the projection transform is applied to a pattern h{x,y) having a variation in only 
one direction with that direction being at an orientation angle >9„ with the horizontal direction, 
then the values of the projection pj^r) is significantly higher when the angle 6 is equal to the 
orientation angle P„ compared to all other values of angle ft Fig. 6 illustrates the projection 
pdr) of the example pattern illustrated in Fig. 4A plotted in the quasi-polar space. It can be 
seen that the values of the projection p^r) are significantly higher when the angle ft is equal to 
the angle fi,,. The projection p^r) also has even or odd symmetry about the offset distance v„. 

Hence, by applying the Radon transfomis on the functions in Equations (17), (18) 
and (19), it can be shown that the Radon transform of such a function is constrained to a line 
where ft = as shown in the following: 

p(r,ft)= RAh{x,y)} = g{r - v,)s{0 - fi) (20) 
(r,ft)= R, {h;^^ {x,y)]^\r - 1./"" 5{0 - /?) (21 ) 

p-(r,ft)= Re{h;^,{x,y)}^sgn{r -v,\r-- r,]'*'" S{e - fi) (22) 
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Having concentrated a two-dimensional function h(x{y) formed from an one- 
dimensional basis function g{r) onto a single line in polar space by the use of Radon 
transform, it is further possible to concentrate the function's energy into a single point (or into 
a small region near a point) by using 1-D correlation detection. If the two-dimensional 
function h{x^) is formed from an one-dimensional basis function g{r) having scale invariant 
properties as herein defined, then the correlation detection will be scale invariant. If the 
watermark includes several basis patterns at different (and in general unknown) angles 
then the basis patterns may be detected by repeating 1-D quasi-radial correlations (in 
coordinate r) for all values of the polar angle 0 . As with other correlations and convolutions, 
a series of 1-D correlations may be efficiently implemented by use of discrete FFT 
algorithms, specifically using 1-D FFTs of 2-D functions along the radial axis. 

In the case of Equation (2 1 ), the 1 -D correlation (symbol ® ) of a 2-D function using 
the dummy radial coordinate r| is as follows: 

{r,0) = |;7 - rf-'V^ {n,0)drj . H^'" <8>p*{r,e) (23) 

—00 

The correlation function c^{i\d) is a highly concentrated function of both r and 9, 
With the preferred parameter p = - 1/2 , the correlation function c^{r,e) has good 
orthogonality properties. For parameter /?«-l the correlation function c*(r,d) resembles a 
near perfect delta spike function 5(r-ro>y(^-/?). A near perfect delta-like result may also 
be advantageously achieved by the use of phase-only-correlation. 

Similarly for the odd form h' (Equation (19)), the 1-D correlation is: 

c-{r.e)= 7sgn(7 - r\i - rf' {ii^e)cin - (sgn(r^ri''^''" )® (24) 

-Oft 

The more general correlation of the pattern h{x,y) expressed in Equation (17) is: 
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c{>;B)= jg{TJ-r)p{tj,0\l,j^g{r)®p{r,0) (25) 

-00 

where the asterisk denotes complex conjugation. 

Step 220 follows where the processor 105 performs a 1-D correlation between the 
projection />6(r) and the basis function in the quasi-radial coordinate r for all possible values 
of the polar angle 0. If the basis function is a complex function, such as that in Equation (1), 
then the complex basis function is preferably correlated with the (real) projection p^r). The 
term "correlation" also includes phase correlation and phase correlation scaled by energy. 
The resulting correlation has a peak at polar coordinates (v«,/?/,) for each embedded pattern. 
The resulting correlation of the 1-D correlation between the projection pc^r) illustrated in Fig. 
6 and the bjisis function of Equation (1) is illustrated in Fig. 7. 

Finally, in step 230 the processor 105 finds the absolute peaks of the correlation. 
The orientation angle J3„ and perpendicular displacement vj, of each embedded pattern is 
directly available from the polar coordinates ( of the peaks. 

Fig. IB shows a flow diagram of a method 400, which is an implementation of the 
method 200 (Fig. lA) of detecting one or more patterns embedded in an image. In this 
implementation, the Radon transform is used as the projective transform, and is implemented 
by resampling in the frequency domain. The correlation with the basis function (step 220, 
Fig. 1 A) is also performed by multiplication in the frequency domain. 

Method 400 starts by the processor 105 performing steps 202 and 205, which are the 
same as those of method 200 described above. 

The processor 105 then computes in step 211 the 2-dimensional Fast Fourier 
transform (FFT) of the image. Preferably, before computing the FFT, the image values 
(intensities) near the image edges are first attenuated by the processor 105 so that the image 
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values fade to zero gradually and smoothly towards the edges of the image. Those skilled in 
the art of image processing are well aware of such preprocessing steps. 

Step 212 follows where the values of the transform, which are in the Cartesian 
frequency space, are resampled by the processor 105 onto the quasi-polar Frequency space. It 
is noted that the angular spacing and the radial scaling need not be constant. A 1 -dimensional 
inverse Fa&t Fourier transform (IFFT) in the radial direction of the remapped transform of the 
image is calculated by the processor 105 in step 213 for all desired angle values. This 
provides the Radon transform of the image, which is equivalent to the output of step 210 of 
method 200 (Fig. lA). 

Method 400 continues in step 214 where a 1 -dimensional FFT in the radial direction 
of the Radon transform is calculated by the processor 105 for all desired angle values. In 
addition, in step 215, the 1 -dimensional FFT of the (one-dimensional) basis function is 
calculated by the processor 105. Alternatively, the basis function may be mathematically 
transformed. The processor 105 next multiplies the transform of the basis function in step 
216 with values of the output of step 214 along constant radial lines for all angle values. The 
IFFT in the radial direction of the output of step 216 is then calculated by the processor 105 in 
step 217. Steps 215 to 217 are the frequency space equivalence of step 220 in method 200 
(Fig. 1 A), and the output of step 21 7 is the same as that of step 220. 

Finally, the processor 105 finds the absolute peaks of the correlation in step 230 in 
the manner described in relation to Fig. 1 A. 

Fig. IC shows a flow diagram of a method 500, which is another implementation of 
the method 200 (Fig, 1 A) of detecting one or more patterns embedded in an image. The steps 
of method 500 are the same as those of method 400 (Fig. IB) except, because steps 213 and 
214 effectively cancel each other out, those steps are not performed in method 500. 
Accordingly, the steps of method 500 have descriptions corresponding to those described 
above. 
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If a pattern is embedded a substantial distance from the centre of an image, or if an 
image is cropped substantially so that the embedded pattern is significantly decentred, then 
the correlation peak may suffer noticeable degradation and asymmetry. Substantial 
degradation can prevent successful detection of an underlying peak. To counteract the 
5 reductions in peak detectability it can be advantageous to use the Hilbert transform detection 
technique. For 1-D scale invariant patterns the Hilbert transform relation is particularly 
simple owing to an unusual self-Hilbert transfonn relationship. In practice the Hilbert 
technique simplifies in this case to correlating the bilateral (i.e two-sided) real embedded 
pattern with two separate complex, and one-sided detection templates. In practice two 

10 complex correlations are carried out with one-sided complex templates that are mirror 
reflections of each other. The square correlation magnitudes are then added together to give 
the overall correlation magnitude squared. The one-side correlations have straightforward 
interpretations in the Fourier domain (as Hilbert transforms) which allows efficient 
implementation by direct use of the Hilbert transforms. 

15 Methods 200, 400 and 500 will detect patterns even if the watermarked image is 

rotated. In order for the detection to detect the patterns after scaling^ anamorphic scaling and 
shear, patterns with scale invariant correlation properties, as described hereinbefore, have to 
be used. 



20 constant factor) for complex functions of infinite extent in the spatial direction. Typically, an 
embedded watermark is realised as a real function of limited extent. The limit in spatial 
extent constrains the scale invariance to be effective within some range, typically from 25% to 
400%. The limit to a real function is more serious in that exact scale invariance, strictly 
speaking, is lost. Consider a real function based on the real part of Equation (1): 



However, true scale invariance is an ideal, and is only achievable (within a complex 



25 




(26) 
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A rescaling by a factor s provides: 



(27) 

= (j|''H''{cos(a In[A:])cos(a ln[s])-sin(a ln[A:])sin(a In[^D} 



new 



ure 



Clearly this is not a simple complex factor as shown in Equation (2). A 
component has been introduced. The new sin(a ln[x]) component is known as the quadrati 
function orthecos(aln[jcj) function. There is an elegant manner for regaining the overall 
scale invariance in the detection methods200, 400 and 500, and that is to use a fully complex 
matching function in the correlation procedure. The general manner can be represented as 
follows: 

<^M= fs'in-rMpM}dn^gir)®^{pir,e]} (28) 

Assuming a general scale invariant spread-spectrum/spread space function g, it can 
be proved that the correlation of a complex matching function with a real embedded function 
gives essentially the same correlation peak as a complex-complex coirelation. but reduced by 
a factor of two. Hence, 



(29) 



There is an insubstantial difference in the detection peak shape, with the background 
difference being orders of magnitude below the peak level. Note that maintenance of the 
line's overall energy (via the Parseval / Plancherel theorem) is important for consistent peak 
values between lines with different values of the orientation angle 0. The energy is 
intrinsically maintained for conventional correlation; however, phase correlation requires 
specific attention to maintain the correct energy for each line. 

As set out above, each pattern included in the embedded watermark has two spatial 
parameters, those being the perpendicular displacement v„ and the orientation angle /}„ . 
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Distortions result in a change in one or both the spatial parameters of the patterns. Fig. 8A 
shows 3 correlation peaks detected in an image having 3 embedded patterns with different 
parameters. Fig. 8B shows the correlation peaks detected from a distorted version of the same 
image. As can be seen when the positions of the peaks in Figs. 8A and SB are compared, the 
spatial parameters of the patterns have changed due to the distortions. 

Rotation, scaling, anamorphic scaling and shear are affine transformations in that 
they preserve straight lines along with the parallelism of lines, but not lengths and angles. 
Such distortions may be represented by a 2x2 matrix operating upon a coordinate vector 
{x,y) to give a coordinate (^',>'') as follows: 



H: :) 



x' = Ax 



(30) 



Rotation may be written in the matrix form of Equation (30) as follows: 



x'\ ( cos P sin p 
y'j \-sinP cosP 



Similarly, scaling may be written as: 



yj 



oY 



0 a 



Ay. 



(31) 



(32) 



Anamorphic scaling along the x and y axes is sometimes called aspect ratio change 
and has the form: 



(33) 



The scaling may be along any two directions, with those directions not necessarily 
being orthogonal. Shear distortion must therefore have the shear direction specified. In the 
case of shear in the x direction, the distortion has the form: 
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y 



0 1 



(34) 



and in the direction: 



(35) 



By using four or more patterns in the watermark, with the patterns having suitable 
5 spatial parameters, the total distortion may be estimated. The total distortion may also be 
decomposed into convenient combinations of the prime distortions, namely rotation, scaling, 
anamorphic scaling and shear. Consider for example a general distortion followed by 
translation, which may be expressed as: 



= 5 



^11 



V^2I ^21 ) 



(36) 



10 wherein 5 is a scaling factor, and {x^^y^) is a translation. When the scaling factor s 

is negative, the distortion is rotation by 180°, Furthermore 



(37) 



indicates that a reflection has occurred. Consider a point (xj/) on the axis of 
symmetry before distortion, the point (j:^) being a distance A„ along the axis of symmetry. 
IS The canonical representation of the point {x;y) is: 



[>' = v,sin>0, +/l,cos^. 



(38) 



Representing the point {xy) as a vector r, and defining a vector p as the unit vector 
perpendicular to the axis of symmetry as follows: 
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r = -;ci + vj 

1 (39) 
[p = cos/?„i + sin/?J 

then the dot product of vectors r and p is equal to the displacement v« as follows: 

r#p = v^ =.vcos>9,, +rsin^„ (40) 

Fig, 4C illustrates the relationship between the vectors r and p, as well as the 
distances v., and \„. Defining cos/?, = c„, and smj3„ = s„, and substituting Equation (38) into 
Equation (36) provides: 



(41) 



jjc' = sa^,{v„c„ - X„s„) + sa„ {u„s„ + A„c„ )+ 
Writing Equation (41) in a different form provides: 

(- «i I J» + «.2C„ )s " (- a,,5„ + iz,,c„ )y 

= -SV^CanC, +«,25„X-a2iJ, +«22«^J- (43) 
•s»'„(a2.c« +ai2*„X-aii*« +«i2<^«) 

Similar to the above, vector f is the vector to point {x\y'), that is point (x^) after 
the distortion, unit vector p is perpendicular to the axis of symmetry of the distorted pattern 
and forming an angle with the horizontal axis, and point ix\y') being a distance X„ along 
the axis of symmetry. The canonical representation of the point {x',y') is: 



x' = v,cosfi„ ~X„sinfi„ 
y = v„sinfi„ +4cos>9„ 



(44) 



Also, the dot product of vectors f and p is equal to the displacement v„ after 
distortion as follows: 
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^ • P = ^, = cos + sin (45) 
Combining Equations (41) and (45) provides: 

^^^p^^^^tA^jLlShl^ (46) 



p _ (gn<^« +^^2^.X-^^2l>y. +^^2gJ-(^^lg. +^^2^»X-^M'y« +^^I2^J 

+«I2^J' +(^'U<^« +^'l2^J' 

V(^II^« +«I2^J' +(a|tC„ +^'l2^«)' 



(47) 



5 The parameters fi^ and i?„ are available from the position of the peak for each 

pattern detected using method 200 (or equivalently method 400 or 500), Equations (46) and 
(47) include 7 unknowns. Therefore, by embedding and detecting at least 4 patterns with 
different parameters, and solving Equation (46) and (47), the distortion parameters r/u, a\2^ 
«2i, «22» xo and>'o can be found. With the distortion parameters an, cfi2, (i2\, ^21, ao and^o 
10 known, the distortions may be reversed to form the undistorted image. 

T>pically, both the perpendicular displacement Vn and the orientation angle fi^ are 

chosen to be widely separated in the embedded patterns to give the best accuracy. Use of 
more than 3 patterns allows noise insensitive redundancy in the arrangement of patterns. 

In genera], more than 3 patterns are required if reflection is also to be detected in 
15 addition to rotation, scaling, anamorphic scaling and shear. Once the distortion parameters 
are known, the original parameters, that is the perpendicular displacement and the 
orientation angle fi^ of each pattern, may be determined. 

Using method 300 for embedding a watermarking and method 200, 400 or 500 to 
detect the embedded watermark, data or a pointer to that data (for example using a URL or 
20 pointer to that URL) may be embedded into an image. The so-called image metadata, that is 
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data about, or referring to some property of that image, is thus bound to the image in that it 
can be retrieved from the image even if the image is distorted. 

The metadata may include a database address, or related data. This allows images to 
be uniquely associated with a database address, even if the image is processed and any 
datafile heading removed through malicious intent or mere format conversion. Similarly, the 
metadata may include copyright infomiation of the image. 

The distortions that the metadata can resist include the projective transforms above, 
but also include the following: 

Printing; 

Photocopying/copying; 

Scanning; 

Colour removal; 

Gamma correction/gamma change; 
JPEG compression/general compression; 
Format conversion (ie BMP to GIF); 
Noise addition and removal; 
Filtering, such as low-pass filtering; and 
Cropping* 

The manner in which the metadata may be encoded into the image includes: 

- Encoding the metadata in an arrangement of embedded patterns, that is, using the 
paramenters /X, and v„ of embedded patterns. The decoding, which follows detection, then 
consists in recognising the arrangement irrespective of projective transformations. 

- Encoding the metadata using an alternative watermarking scheme. Additional 
patterns are embedded using method 300, which are then detected using the method 200, 400 
or 500. Using the additional patterns, the image distortion is then determined, followed by 
reversing the distortion. The watermarks used to encode the metadata may then be detected 
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and decoded in the usual manner Alternatively, the transformation, once determined, may be 
used to define a mapped detection of the alternative watermark, followed by a mapped 
decoding, if required. This approach avoids direct image inversion and the associated 
computation of image warping/interpolation. 
5 In some situations it is not necessary that the embedded patterns be imperceptible to 

the human eye. Nor may it be necessary that the object undergoing embedding be an image. 
For example a planar object, such as a silicon wafer, may have a number of patterns 
embedded, using engraving or etching, in its surface. The patterns may then be and utilised to 
determine the surface's position and orientation. 

10 Fig. 1 1 shows a schematic block diagram of a system 600 that may be used to 

implement the embedding and detection of watermarks in an image. The system 600 includes 
a first apparatus 610 for receiving a first image 601, and for embedding a watermark into the 
image 601 to produce a second image 602. 

The first apparatus 610 may be the general-purpose computer system 100 described 

15 in relation to Fig. 3 wherein the first image 601 is captured using imaging device 122, or 
received firom the network 120. In this case, the second image 602 may be printed on printer 
1 15 or directly communicated to another apparatus, such as a second apparatus 620, over the 
network 120. 

In alternative implementation the method 300 is performed directly in a scanner or a 
20 camera, where the first image 601 is captured using an imaging device of the scanner or 
camera, and the second image 602 is in electronic form. 

The system 600 further includes a second apparatus 620 for receiving a third image 
621 that may contain an embedded watermarks, and for detecting the watermark, if such a 
watermark is present. The third image 621 may be the same as the second image 602, or may 
25 be a distorted version of the second image. 
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The second apparatus 610 may also be the general-purpose computer system 100 
described in relation to Fig. 3 wherein the third image 601 is captured using imaging device 
122, or received from the network 120. 

The second apparatus 620 is operative to then detect the watermark, if a watermark 
exists, and report such to a user. The second apparatus 620 may additionally decode any 
information that may be embedded in the placement of the respective patterns making up the 
watermark. In a case where the second apparatus 620 expect to find an image with embedded 
patterns having predetermined parameters, the second apparatus 620 may detect that the third 
image 621 is a distorted version of the second image, and such distortions may be 
automatically reversed to produce a fourth image 622. 

In an alternative implementations one of methods 200, 400 or 500 is performed 
directly in a scanner, where the third image 62 1 is captured using the scanner. The scanner 
may then alert the user, may automatically reverse distortions to produce the fourth image 
622, or decode information encoded in the parameters of the embedded patterns. 

The foregoing describes only some embodiments of the present invention, and 
modifications and/or changes can be made thereto without departing from the scope and spirit 
of the invention, the embodiments being illustrative and not restrictive. 

In the context of this specification, the word "comprising" means 'Including 
principally but not necessarily solely" or "having" or "including", and not "consisting only 
of. Variations of the word "comprising", such as "comprise" and "comprises" have 
correspondingly varied meanings. 
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The claims deflning the invention are as follows: 



1. A method of detecting one or more patterns embedded in an image, each pattern 
having been formed substantially from a one-dimensional basis function, said method 

. 5 comprising the steps of: 

(a) calculating a projective transform of said image; 

(b) calculating a 1-D correlation between the projective transform and said 
basis function for a selection of angles; and 

(c) finding one or more peaks of said correlation, wherein the position of each 
10 of said peaks provides spatial parameters of one of said one or more embedded patterns. 

2. - A method of determining transformations applied to an image, said image having 
at least three patterns embedded therein, each pattern having been formed substantially from a 
one-dimensional basis function being substantially scale invariant as herein defined, said 

15 method comprising the steps of: 

(a) calculating a projective transform of said image; 

(b) calculating a 1-D correlation between the projective transform and said 
basis function for a selection of angles; 

(c) finding peaks of said correlation, wherein the positions of said peaks 
20 provide spaiial parameters of said embedded patlems; and 

(d) determining from said spatial parameters said transformations. 



3. A method as claimed in claim 1 or 2, wherein step (a) comprises the sub-steps of: 

(al) transforming said image to the frequency domain using a two-dimensional 

25 transform; 

(a2) resampling the transformed image onto a quasi-polar map; and 
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(a3) inverse transforming the resampled image to the space domain using a one- 
dimensional transform. 

4- A method as claimed in any one of claims 1 to 3, wherein step (b) comprises the 

sub-steps of: 

(bl) transforming said projective transform to the frequency domain using a one- 
dimensional transform; 

(b2) transforming said basis function to the frequency domain using a one- 
dimensional transform; 

(b3) multiplying the transformed projective transform with the transformed basis 
function along radial lines to create a result; and 

(b4) inverse transforming said result to the space domain using a one- 
dimensional transform. 

5. A method of detecting one or more patterns embedded in an image, each pattern 

having been formed substantially from a one-dimensional basis function, said method 
comprising the steps of: 

transforming said image to the frequency domain using a two-dimensional 

transform; 

resampling the transformed image onto a quasi-polar map; 

transforming said basis function to the frequency domain using a one-dimensional 

transform; 

multiplying the resampled transformed image with the transform of said basis 
function along radial lines to create a first result; 

inverse transforming said first result to the space domain using a one-dimensional 
transform to create a second result; and 
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finding one or more peaks of said second result, wherein the position of each of 
said peaks provides spatial parameters of one of said one or more embedded patterns. 

6. A method of determining transfomiations applied to a first image, said first image 

being a transformed copy of a second image, said second image having at least three patterns 
embedded therein, each pattern having predefined parameters and having been formed 
substantially from a one-dimensional basis function being substantially scale invariant as 
herein defined, said method comprising the steps of: 

calculating a projective transform of said first image; 

calculating a 1-D correlation between the projective transform and said basis 
function for a selection of angles; 

finding peaks of said correlation, each peak corresponding with one of said 
embedded patterns; and 

determining said transformations from the positions of said peaks and said 
parameters of said patterns embedded into said second image. 

7. A method of embedding a watermark into an image, said method comprising the steps 
of: 

maintaining a basis function, said basis function being a scale invariant one- 
dimensional function excluding the function /(/-)- cos(a log|r| + c) wherein a and c are 
constants; 

forming one or more patterns from said basis function, each pattern having variation in 
one direction; and 

adding said basis pattem(s) to said image. 
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8. A method as claimed in claim 7, wherein said scale invariant one-dimensional 
function is a function wherein the profile of the phase of said function is substantially 
invariant after scaling within a predefined range of values. 

5 9. A method as claimed in claim 8, wherein said scale invariant one-dimensional 
function is selected from the group comprising: 

/, ('•) = [sgn(r)]* exp{£(ar-^ + e)); 

10 /3O')- [sgn(r exp^/— A-^ j; wherein y, £, p, M, , a, k and ^ are constants 

and is an arbitrary function. 

10. A method as claimed in any one of claims 7 to 9, wherein said basis function is a 
complex function and said one or more basis patterns are real valued. 

15 

11. A method as claimed in any one of claims 7 to 10, wherein each pattern is added with 
a respective displacement from the image centre and a respective orientation. 

12. An apparatus for detecting one or more patterns embedded in an image, each 
20 pattern having been formed substantially from a one-dimensional basis function, said 

apparatus comprising: 

means for calculating a projective transform of said image; 
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means for calculating a 1-D correlation between the projective transform and said 
basis function for a selection of angles; and 

means for finding one or more peaks of said correlation, wherein the position of 
each of said peaks provides spatial parameters of one of said one or more embedded patterns. 

5 

13. An apparatus for determining transformations applied to an image, said image 
having at least three patterns embedded therein, each pattern having been formed substantially 
from a one-dimensional basis function being substantially scale invariant as herein defined, 
said apparatus comprising: 

10 means for calculating a projective transform of said image; 

means for calculating a 1-D correlation between the projective transform and said 
basis function for a selection of angles; 

means for finding peaks of said correlation, wherein the position of said peaks 
provide spatial parameters of said embedded paiiems; and 
15 means for determining from said spatial parameters said transformations. 

14. An apparatus as claimed in claim 12 or 13, wherein said means for calculating a 
projective transform of said image comprises: 

means for transforming said image to the frequency domain using a two- 
20 dimensional transform; 

means for resampling the transformed image onto a quasi-polar map; and 
means for inverse transforming the resampled image to the space domain using a 
one-dimensional transform. 

25 15. An apparatus as claimed in any one of claims 12 to 14, wherein said means for 

calculating a 1-D correlation comprises: 
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means for transforming said projective transform to the frequency domain using a 
one-dimensional transfonn; 

means for transforming said basis function to the frequency domain using a one- 
dimensional transform; 

means for multiplying the transformed projective transform with the transformed 
basis function along radial lines to create a result; and 

means for inverse transforming said result tQ the space domain using a one- 
dimensional transform. 

16. An apparatus for detecting one or more patterns embedded in an image, each 
pattern ha^'ing been formed substantially from a one-dimensional basis function, said 
apparatus comprising: 

means for transforming said image to the frequency domain using a two- 
dimensional transform; 

means for resampling the transfoniied image onto a quasi-polar map; 

means for transforming said basis function to the frequency domain using a one- 
dimensional transform; 

means for multiplying the resampled transformed image with the transform of 
said basis fimction along radial lines to create a first result; 

means for inverse transforming said first result to the space domain using a one- 
dimensional transform to create a second result; and 

means for finding a peak of said second result, wherein the position of each of 
said peaks provides spatial parameters of one of said one or more embedded patterns. 

17. An apparatus for determining transformations applied to a first image, said first 
image being a transformed copy of a second image, said second image having at least three 
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pattems embedded therein, each pattern having predefined parameters and having been 
fonned substantially from a one-dimensional basis function being subslantialiy scale invariant 
as herein defmed, said apparatus comprising: 

means for calculating a projective transform of said first image; 

means for calculating a 1-D correlation between the projective transform and said 
basis function for a selection of angles; 

means for finding peaks of said correlation, each peak corresponding with one of 
said embedded patterns; and 

means for determining said transformations from the positions of said peaks and 
said parameters of said patterns embedded into said second image. 

18. An apparatus for embedding a watemiark into an image, said apparatus comprising: 
means for maintaining a basis function, said basis function being a scale invariant one- 
dimensional function excluding the function /(r) = cos(a log|r| + c) wherein a and c are 
constants; 

means for forming one or more patterns from said basis function, each pattern having 
variation in one direction; and 

means for adding said basis pattem(s) to said image. 

19. An apparatus as claimed in claim 18, wherein said scale invariant one-dimensional 
function is a function wherein the profile of the phase of said function is substantially 
invariant afler scaling within a predefined range of values. 

20. An apparatus as claimed in claim 19, wherein said scale invariant one-dimensional 
function is selected from the group comprising: 
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X{f) = >-''exp(i-r''\ 
V r ) 

/, (r) = [sgn(/-)]* expiiica -" + s)); 
Airh 2^„/o j//|V); 



\ r J 

and fo is an arbitrary function. 



; wherein /?, N\, N2, a, * and /u are constants 



21 . An apparatus as claimed in any one of claims 18 to 20, wherein said basis fiinction is a 
complex function and said one or more basis patterns are real valued. 

22. An apparatus as claimed in any one of claims 18 to 21, wherein each pattern is added 
with a respective displacement from the image centre and a respective orientation. 

23. A program stored in a memory for detecting one or more patterns embedded in an 
image, each pattern having been formed substantially from a one-dimensional basis function, 
said program comprising: 

code for calculating a projective transform of said image; 

code for calculating a 1-D correlation between the projective transform and said 
basis function for a selection of angles; and 

code for finding one or more peaks of said correlation, wherein the position of 
each of said peaks provides spatial parameters of one of said one or more embedded patterns. 

24. A program stored in a memory for determining transformations applied to an 
image, said image having at least three patterns embedded therein, each pattern having been 
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formed substantially from a one-dimensional basis function being substantially scale invariant 
as herein defined, said program comprising: 

code for calculating a projective transform of said image; 

code for calculating a 1-D correlation between the projective transform and said 
5 basis function for a selection of angles; 

code for finding peaks of said correlation, wherein the positions of said peaks 
provide spatial parameters of said embedded patterns; and 

code for determining from said spatial parameters said transformations. 

10 25. A program stored in a memory for detecting one or more patterns embedded in an 

image, each pattern having been formed substantially from a one-dimensional basis function, 
said program comprising: 

code for transforming said image to the frequency domain using a two- 
dimensional transform; 
15 code for resampling the transformed image onto a quasi-polar map; 

code for transforming said basis function to the frequency domain using a one- 
dimensional transform; 

code for multiplying the resampled transformed image with the transform of said 
basis function along radial lines to create a first result; 
20 code for inverse transforming said first result to the space domain using a one- 

dimensional transform to create a second result; and 

code for finding one or more peaks of said second result, wherein the position of 
each of said peaks provides spatial parameters of one of said one or more embedded patterns. 



25 26. A program stored in a memory for determining transformations applied to a first 

image, said first image being a transformed copy of a second image, said second image 
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having at least three patterns embedded therein, each pattern having predefined parameters 
and having been formed substantially from a one-dimensional basis function being 
substantially scale invariant as herein defined, said program comprising: 

code for calculating a projective transform of said first image; 
5 code for calculating a 1-D correlation between the projective transform and said 

basis function for a selection of angles; 

code for finding peaks of said conrclalion, each peak corresponding with one of 
said embedded patterns; and 

code for determining said transformations from the positions of said peaks and 
10 said parameters of said patterns embedded into said second image. 

27. A program stored in a memory for embedding a watermark into an image, said 
program comprising: 

code for maintaining a basis function, said basis function being a scale invariant one- 
15 dimensional function excluding the function /(r) = cos(a log|r| -h c) wherein a and c are 

constants; 

code for forming one or more patterns from said basis function, each pattern having 
variation in one direction; and 

code for adding said basis pattem(s) to said image. 

20 

28. A program as claimed in claim 27, wherein said scale invariant one-dimensional 
function is a function wherein the profile of the phase of said function is substantially 
invariant after scaling within a predefined range of values, 

25 29. A method substantially as herein described with reference to the accompanying 

drawings. 
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30. An apparatus substantially as herein described with 

accompanying drawings. 

DATED this Third Day of October 2002 
Canon Kabushiki Kaisha 
Patent Attorneys for the Applicant 
SPRUSON&FERGUSON 
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